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ABSTRACT 

The presence of substructure in galaxy groups and clusters is believed to be a sign 
of recent galaxy accretion and can be used not only to probe the assembly history of 
these structures, but also the evolution of their member galaxies. Using the Dressler- 
Shectman (DS) Test, we study substructure in a sample of intermediate redshift 
(z ~ 0.4) galaxy groups from the Group Environment and Evolution Collaboration 
(GEEC) group catalog. We find that 4 of the 15 rich GEEC groups, with an average 
velocity dispersion of ^525 km s _1 , are identified as having significant substructure. 
The identified regions of localized substructure lie on the group outskirts and in some 
cases appear to be infalling. In a comparison of galaxy properties for the members of 
groups with and without substructure, we find that the groups with substructure have 
a significantly higher fraction of blue and star-forming galaxies and a parent colour 
distribution that resembles that of the field population rather than the overall group 
population. In addition, we observe correlations between the detection of substructure 
and other dynamical measures, such as velocity distributions and velocity dispersion 
profiles. Based on this analysis, we conclude that some galaxy groups contain sig- 
nificant substructure and that these groups have properties and galaxy populations 
that differ from groups with no detected substructure. These results indicate that the 
substructure galaxies, which lie preferentially on the group outskirts and could be in- 
falling, do not exhibit signs of environmental effects, since little or no star-formation 
quenching is observed in these systems. 

Key words: galaxies:statistics-galaxies:kincmatics and dynamics-galaxies:formation 



1 INTRODUCTION 

The current theory of structure formation in the Universe is 
based on the standard A cold dark matter (ACDM) cosmo- 
logical model, in which objects grow hierarchically from the 
initial matter density perturbations through mergers and 
accretion fe.g..|Press fc Schechter|[l97i ; lLacev fe Colli 19931 ; 
ISpringel et al.| [2005h In order to test the theory of hierarchi- 
cal structure formation, one must investigate the assembly 
history of the large structures in the Universe, namely galaxy 
groups and clusters. 

A natural consequence of a hierarchical Universe is the 
existence of substructure within larger systems. Tradition- 
ally, substructure has been defined as a kinematically dis- 
tinct sub-halo within a larger parent halo. A broader, more 
observationally motivated definition, and one that we will 



assume here, also includes separate haloes that are either 
merging to form a larger halo; gravitationally bound and 
infalling onto a pre-existing halo; or in the nearby large- 
scale-structure, but not necessarily bound or infalling. In 
general, accreting structure does not have the same kine- 
matic properties as the host, but whether or not substruc- 
ture can be observed depends on how long it remains in- 
tact after infall. Early iV-body simulations suggested that 
the assimilation of substructure into the h ost was rapid, 
providing no detectable lo ng-lived features l|Katz fc White! 
1 19931 ; IS ummers et alJll995h . However, later work has shown 
that the lack of observable su bstructure in these simula- 
tions was due to poor resolution (I Moore et al Il996lh Indeed, 



high resolution A^-body s imulations (e.g., iDiemand et al.l 
2008; ISpringel ct al. 2008), demonstrate that several hun- 
dred thousand sub-haloes can exist in a Milky Way sized 
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dark matter halo at a redshift of zero. Using semi-analytic 
models to study the substructure within ind ividual galaxy to 
cluster-sized haloes, iTavlor fc Babul (|2004l ) showed that ac- 
creting systems could survive within the host halo for many 
orbits, depending on the orbital parameters of the substruc- 
ture upon infall. 

These theoretical results suggest that substructure 
should be a detectable quantity and numerical dark 
matter simulations of galaxy groups and clusters in a 
ACDM Universe predict that approximate ly 30 per cent of 
all sy stems should contain substructure jKnebe fc M uller 

(e.g., iBeers et al 



2000). Studies of individual clu s ters , , 
1982 : iDressler fc Shectmanl Il988l; I West fc Bothunl 
Bird! Il994al; IColless fc Dunnl 1 19961 : burgett fc et ail 



1990 



2004 



Bohringer et al. 2010h indicate that a large fraction of sys- 



tems show evidence of significant sub-clustering. The pre- 
dicted theoretical value of 30 per cent is in agreement 
with some substructure studies of groups and clu sters (e.g., 
IZabludoff fc Mulchaevl[T998al : ISolanes et ailll999l i but sev- 
eral other result s have demon strated a much higher fraction 
of substructure; iBir 3 (|l994bh observe d that 44 per cen fl of 
their sample contained substruc ture, [D ressier & Sh ectmanl 
(|l988h observed 53 per cent and iRamella et all (|2007T ) find 
an extremely high substructure fraction of 73 per cent. Al- 
though the precise fraction of groups and clusters with sub- 
structure may still be a source of debate, the observed pres- 
ence of any substructure strongly suggests that these sys- 
tems grow in a hierarchical manner through the accretion of 
galaxies and smaller groups of galaxies. 

Though galaxy clusters are amongst the largest struc- 
tures in the local Universe, they do not represent the 
most common host environment for galaxies. Galaxy groups, 
which contain a few to tens of member galaxies, are the host 
environmen t of more than half of the present-da y galaxy 
population dGeller fc Huchrall 19831 : lEke et al.ll2005h . Despite 
the importance of groups in the build up of large galaxy clus- 
ters, there have been few studies on the assembly history of 
groups themselves. 

One method of probing group assembly histories is to 
look at the amount of substructure within these systems. 
The presence of such structure would indicate that th e group 
has recently accreted galaxies (lLacev fc Cold 1 19931 ). Stud- 
ies have be en carried out for gala x y grou ps in the local 
Universe by IZabludoff fc Mulchaevl ( 1998a), who observed 
significant (> 99 per cent confidence level) substructure in 
two of their six local groups. An interesting result of their 
analysis showed that the substructure was located on the 
outskirts of the systems, that is ~0.3-0.4 h" 1 Mpc, where 
h = ffo/(100 km s~ x Mpc -1 ), from the core of the g roup. 



Based on these results, IZabludoff fc Mulchaevl l|l998al ) con- 
cluded that like rich galaxy clusters, some galaxy groups 
assembled in a hierarchical manner through the accretion of 
smaller structures fr om the field . A m ore recent study of lo- 
cal galaxy groups bv lFirth et al.l (|2006l ) found similar results, 
with roughly half of their sample s howing significant sub- 
structu re. A lthough, the findin gs of IZabludoff fc Mulchaevl 
jl998al ) and lFirth et"ail (2006) provided an important first 
step in unveiling the assembly history of groups, their anal- 



1 11/25 galaxies in their sample have significant (> 95 per cent) 
substructure based on results of the Dressler-Shectman statistics. 



ysis was based on a small sample of very nearby systems. In 
order to gain a better understanding of the role of groups in 
the growth of structure, we must search for substructure in 
a larger sample of galaxy groups with highly complete spec- 
troscopy that cover a wide range in redshift. Such studies 
have only become possible with recent deep spectroscopic 
surveys that have produced large gro up catalogues, such as 
the Sloan Digital Sky Survey (SDSS) (jBerlind fc et al.ll2006l : 
lYang et al.l 120071 ) , the Group Environment and Evolution 
Colla boration (GEEC) (|Wilman et al.ll2005l : ICarlberg et all 
l200ll). and the high er redshift extension of GEEC (GEEC2) 
I Balo gh et al.ll201ll ) optical group catalogs. 

In addition to the role of the group environment in the 
growth of large scale structure in the Universe, studies of 
substructure within groups allows us to probe galaxy evolu- 
tion. Since groups have fewer members than galaxy clusters, 
any substructure present will have a stronger effect on the 
observed group properties (i.e. colours, blue or active frac- 
tions). If substructure traces accreting galaxies, one might 
expect a correlation between observed galaxy properties and 
substructure. Possible correlations could exist between sub- 
structure and the colours of galaxies in groups, or with sub- 
structure and star formation rates. One of the main goals of 
this paper is to search for such correlations. 

In this paper we search for substructure in a sample of 
intermediate redshift galaxy groups from the GEEC Group 
Catalog. In Section 2, we discuss the sensitivity and reli- 
ability of the Dressler-Shectman (DS) Test for group-sized 
systems using Monte Carlo simulations. In Section 3, we ap- 
ply the DS Test to the GEEC groups and discuss the results 
of our analysis. In Section 4, we discuss the relationship be- 
tween the presence of substructure and other indicators of 
dynamical state, such as the shape of the group velocity 
distribution. In Section 5, we look for correlations between 
substructure and the properties of members galaxies, such 
as colour and star formation rate, and discuss the possi- 
ble implications of our findings. In Section 6 we present our 
conclusions. Additionally, we include an Appendix, which 
provides detailed results of our Monte Carlo Simulations. 

Throughout this paper we assume a ACDM cosmology 
with Q.m = 0.3, fl A = 0.7, and H = 75 km s -1 Mpc -1 . 



2 DETECTING SUBSTRUCTURE IN GROUPS 

Numerous tests for substructure have been developed and 
carrie d out for cluster-sized systems (see iPinknev et al.l 
1996, for a thorough review). In a comparison of five 3- 
dimensional (3- D) tests, which use b oth velocity and spa- 
tial information, IPinknev et al.l (| 19961) determined that the 
Dressler-Shectman (DS) Test (|Dressler fc Shectmanl Il988l ) 
was the most sensitive test for substructure, for systems with 
as few as 30 members. In the following section we look at 
the DS Test, and determine its reliability and robustness for 
smaller group-sized systems. 



2.1 The Dressler-Shectman (DS) Test 

Substructure manifests itself as detectable deviations in the 
spatial and/or velocity structure of a system. The aim of 
the DS Test is to compute local mean velocity and velocity 
dispersion values, for each individual galaxy and its nearest 
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neighbours, and then compare these to the global group va l- 
ues. Following the notation of lDressler fc Shectmanl (|l988h . 
we define (p, a) as the mean velocity and velocity disper- 
sion of the entire group, which is assumed to have n mem b er s 
galaxies. Then for each galaxy i in the group, we select it 
plus a number of its nearest neighbours, 7V nn , and compute 
their mean velocity v^^i and velocity dispersion crJ ocal . From 
these we compute the individual galaxy deviations, Si, given 
as 



6 = 



jV„„ + l 



^local -V] + 



°"local 



(1) 



where 1 < i < n mem bers (i.e. Si is computed for each galaxy 
in the system). iDressler fc Shectmanl (|l988T ) originally de- 
veloped the test for cluster-sized systems and the number of 
nearest neighbours used to compute the Si values was rel- 
atively high (i.e. N nn = 11). Since substructure in galaxy 
groups is likely to have fewer than 11 constituent galax- 
ies, we take iV nn to be ^/n mcm b C rs following the methodol- 
ogy of authors who have previously applied the DS Test 
to group-size systems (fsilvermanll 19861 ; IPinknev et al.|[l99^ ; 
IZabludoff fc Mulchaevlll998bl ). This ensures that large kine- 
matic deviations of a few neighbouring galaxies do not 
get diluted by adding too many unassociated galaxies, and 
thereby lowering the computed z7i oca i and o"i oca i values in 
Equation [T] 

The statistic used in the DS Test is the A-value, given 

by 



(2) 



A system is then co nsidered to have substruc ture if 
A/n me mbers > 1.0 (|Dressler fc Shectmanl Il988h . This 
method of using a threshold value to find substructure is 
referred to as the critical values method. 

An alternative method of identifying substructure with 
the DS Test is to use probabilities (P-values) rather than 
critical values. The P-values for the DS Test are computed 
by comparing the observed A-value to 'shuffled' A-values, 
which are computed by randomly shuffling the observed ve- 
locities and re-assigning these values to the member posi- 
tions, a procedure called 'MC shuffling'. The P-values are 
given by 



P — (Ashufflcd > A bs, 



I /n s 



(3) 



where A s huffl e d and A b aer ved are both computed using Equa- 
tion[2]and n s h u ffl c is the number of 'MC shuffles', and there- 
fore the number of A s h u fn e d-values, used to compute the 
probability . One can see from Equation[3]that systems with 
significant substructure will have low P-values, since it is un- 
likely to obtain the observed A-value randomly (Equation 

HJ. 

2.2 Monte Carlo Simulations 



Although IPinknev et al.l (|l996T ) have carried out an exten- 
sive investigation of the DS Test, their analysis was per- 
formed on systems with a minimum of 30 members, and 
the majority of our intermediate redshift groups have fewer 
member galaxies (see Section 3.1). Therefore, we perform 



our own Monte Carlo simulations in order to check the re- 
liability, sensitivity and robustness of the DS Test, using 
both the critical value and probability (P-value) methods, 
for group-sized systems. It should be noted that we specifi- 
cally model our host mock groups after the observed GEEC 
group sample. 



2.2.1 Generating the Mock Groups 

Mock galaxy groups, both with and without substructure, 
are generated using Monte Carlo methods to assign member 
galaxy positions. These mock groups are then used to com- 
pute the false negative and positive rates of the DS Test, and 
also determine the sensitivity of the test against a variety of 
input parameters. 

The radial positions for the members of the mock groups 
are randomly drawn from fits to the projected group-centric 

radial distributions of the galaxies observed in the GEEC 

i 1 

group catalog (jWilman et al.H 2005). The groups are divided 

into four bins based on the number of members in each group 

(Fig. [I}, and fits to each bin are used to populate the mock 

groups. Since the GEEC groups span such a wide range of 

masses and group membership, we elect not to use a single 

fit to the radial distribution of all the group galaxies in our 

sample. Instead, we divide our sample into bins of group 

membership (5 < n mcmbc r S < 10, 10 < n mcmb cr S < 15, 15 < 

^members < 20 and n mem bcrs > 20) and fit each distribution 

separately. It should be noted that we bin our sample by 

group membership, rather than mass or velocity dispersion, 

as we aim to study the false positive and negative rates of the 

DS Test as a function of n me mbcrs- However, the results are 

similar if we bin by a rather than n mcm bcrs- The histograms 

and fits to the radial distributions of the binned groups are 

shown in Fig. [T] The general form of the radial distribution 

of member galaxies is given by 



N oc exp (-XR) . 



(4) 



where R is the radial position of the given galaxy and A = 
2.98, 1.31, 0.902 and 0.606 Mpc" 1 for the 5 < n mcmbc rs < 10, 

10 < ^members ^ 15, 15 < ^members ^ 20 and ^members ^ 20 

bins, respectively. From these results it is clear that groups 
with fewer members have steeper radial distributions and 
smaller maximum group centric radii (Fig. [TJ. 

The redshifts of the member galaxies in the mock groups 
are randomly drawn from Gaussian distributions. The dis- 
persions of the input Gaussians are taken to be the average 
velocity dispersion of the aforementioned group membership 
bins, which are: 300, 350, 400 and 550 km s" 1 . These disper- 
sion values are chosen for the mocks groups in order to best 
mimic the systems in our GEEC sample. We also generate 
mock groups with the same input host velocity dispersion 
("host) f° r all values of n mem bers and find no significant dif- 
ference in our results (see Appendix for further detail). 

As will be shown in Section 2.2.2, in order to determine 
the false negative rates for the DS Test, we must include sub- 
structure within the host group. The positions and redshifts 
of the substructure galaxies are drawn from a Gaussian dis- 
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Figure 1. Histograms of the radial distribution of the galaxies 
in the GEEC Group Catalog stacked in bins of ^members- 

The 

dashed line corresponds to an exponential fit to the distribution, 
which we use to generate the radial positions of galaxies in our 
mock groups. 

tributiorjf]. The free parameters in our mock groups with 
substructure are: the velocity dispersion of the substructure 
(o"redshift), spatial dispersion of the substructure (a pos ition), 
location of the substructure, in both position (e p0 sition) and 
velocity (e re dshift) space, the number of galaxies in the sub- 
structure (n su b), the number of galaxies in the host group 
(^members), and velocity dispersion of the host group (a host ). 
We briefly analyze the effects of each of these parameters on 
the robustness of the DS Test in the following sections, and 
give a more detailed discussion in the Appendix. 

2.2.2 Testing for False Positives 

Using the mock groups with no input substructure described 
in Section 2.2.1, we determine the false positive rates, or type 

1 errors, which for the DS Test occur when substructure is 
identified when none exists. The only parameter we vary for 
these systems is the number of member galaxies (rimcmbcrs)- 

We first compute the false positive rates using the crit- 
ical values method. For each value of n mom b C rs we com- 
pute the A-statistic, (Equations [T] and [2]| and then clas- 
sify a group as having substructure if A/n me mbors > 1.0. 
With this criterion, we find that for all relevant values of 
^members the false positive rates are extremely high. For 
example for n mem bers = 20, the false positive rate is ~81 
per cent and even for larger systems, n me mbcrs = 100, we 
find that the rate is equal l y high . A similar result was ob- 
served bv lKnebe fc Mulled (|2000t ) in a study of substructure 
in numerically simulated galaxy clusters. They found that 
for haloes with no substructure the A/n me mbors- values were 

2 Tests with non-Gaussian input substructure have also been car- 
ried out and for all reasonable distributions the results do not 
differ significantly from the results presented here. 



often greater than 1.0, with values peaking closer to 1.4 for 
larger (i.e. n mcm b C rs ~80 -100) clusters. Although it is possi- 
ble to better identify substructure in richer systems, such as 
clusters, with a higher value of A /remembe rs ~1.4 - 1.6 (e.g., 
iGirardi et af1l2005l ; iKnebe fc Mulledl200d) . we find that for 
group-sized systems this methodology could not simultane- 
ously produce low false positive and false negative rates. 

Alternatively, when P-values are used to identify sub- 
structure, we find that the false positive rates are much 
lower than with the critical values method and also re- 
markably stable for a wide range of group members (i.e. 
5 < n mcm bcrs < 50). From our mock groups we find that for 
significance levels of 0.01, 0.05 and 0.10, the false positive 
rates are 5, 10 and 15 per cent for all values of n me mbers 
tested. In other words, the DS Test will falsely identify 
roughly 5 per cent more substructure than the desired sig- 
nificance level. Although these rates are higher than the ex- 
pected values for a given significance level, they are still 
substantially lower than the rates obtained using the criti- 
cal values method. Also, as long as one chooses significance 
levels of 0.01 or 0.05, the rate of false identifications is ac- 
ceptably low. Based on these results, we rule out the critical 
values method and only perform the following analysis using 
the probabilities, P-values, method. 



2.2.3 Testing for False Negatives 

Having determined the rate of false identifications obtained 
with the DS Test, we now test for the rate of false negatives, 
or type II errors. For this statistic the false negative rate 
measures how often the DS Test fails to detect substructure 
when indeed it exists. The purpose of conducting these false 
negative tests is twofold: first, they help determine the relia- 
bility of the statistic and second, by varying only one of the 
free input parameters (listed in Section 2.2.1) at a time, you 
can determine how sensitive the test is to each parameter. 
The latter places quantitative constraints on the maximum 
size and location of substructure that can be identified. 

Before we begin looking at the input parameters, we 
first determine how reliable the test is at identifying 'ob- 
vious' substructure, that is to say galaxies that are tightly 
correlated and located far from the group centre, both in 
projection on the sky and along the line-of-sight (see Table 
[T] for input parameter values). For these mock groups, we 
find that substructure is correctly identified in almost all 
cases, using the P- value method and a significance level of 
0.01, producing false negative rates of or 1 per cent. 

We then investigate the individual free parameters in 
more detail to determine how each alters the false negative 
rate. We briefly discuss the main results here and leave the 
detailed quantitative analysis, as well as full tables of false 
negative rates, for the Appendix. 

Of the five free parameters, we find that the DS Test 
is most sensitive to the number of galaxies in the substruc- 
ture (n su b) and the location of the peak of the substruc- 
ture's velocity distribution with respect to the peak of the 
host's (e r ed shift). Previous studies of sub stru cture in galaxy 
cluste rs bv lDressler fc Shectmanl l|l988h and iPinknev et all 
(1996) showed that the DS Test was unable to find sub- 
structure that was 'superimposed' with the highest density 
regions of the h osts, s ince the galaxies were spatially mixed. 
IPinknev et"all i|l996l ) even stated that in these cases, any 
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Table 1. Input parameter values for the 'base' mock groups described in Section 2.2.3; i.e. Monte Carlo groups with a zero false negative 
rate. 
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km s" 1 
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0.01 
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15 


5 


400 


0.01 


0.5 


100 


1300 


20 


5 


550 


0.01 


0.5 


100 


1300 


50 


10 


550 


0.01 


0.5 


100 


1300 



type of 3-D test would not be able to accurately detect sub- 
structure. An important distinction, not made by these au- 
thors, is that two forms of 'superposition' can occur. Sub- 
structure can be superimposed with the host group either in 
projected angular position or in redshift space. Our analy- 
sis shows that the DS Test is significantly more sensitive to 
ErcdsMft than to Eposition- The test can usually identify sub- 
structure with superpositions as projected on the sky, but 
has a very difficult time with those along the line-of-sight. 
From Equation [T] it is clear that only collections of neigh- 
bouring galaxies with a different local mean velocity and 
velocity dispersion will produce high Si values, no matter 
their angular position. 

In addition we also find that the level of sensitivity of 
the DS Test to e re dshift) the distance between the substruc- 
ture's and host's peak velocity distribution, is dependent on 
the number of members in the host group. We find that for 
small groups (i.e. n mcm bers < 20), the input substructure 
needs to be further than 2<Thost from the group centre in 
order to be detected. On the other hand, the DS Test can 
identify substructure that is roughly 2<Thost away from the 
peak of the host's Gaussian velocity distribution in groups 
with more than 20 members. For even larger systems, that 
is clusters with n mcm b crs > 50, the input substructure can 
be located within lohost and still be identified. 

Another result from this analysis is that the DS Test 
is sensitive to the number of galaxies that are part of the 
substructure (n su b). Our simulations show that no matter 
the membership of the host group, the test cannot identify 
substructure with fewer then four members. Also, the mini- 
mum required number of members within the substructure 
increases with n m0 mbcrs- This is due to the fact that we set 
N„„ in Equation [T] to -^/n mom bors- Therefore, more members 
in the host group means that the velocity information of 
more 'neighbours' will be used to calculate V\ ocal and a\ ocal . 
Thus, if there are too few galaxies in the substructure, their 
kinematic deviations can be washed out by other neighbour- 
ing host galaxies. 

We find that for rich galaxy groups, with n mem bors > 
20, the DS Test can reliably identify true substruc- 
ture, as has been shown by se veral other authors 
(e.g.. iDressler fc Shectmanl Il988l ; iPinknev et all Il99fj ; 
iKnebe fc Mullerl |2000| ).~ For systems with fewer member 
galaxies, such as poor groups with n mcm b crs < 20, the false 
negative rates are very low (< 5 per cent) only for tightly 
correlated substructure galaxies with large kinematic devia- 
tions. In these poor groups, substructure that is more loosely 
associated with a velocity peak close to that of its host group 
is not as easily detected by the test. 



Taking into account the results of both the false neg- 
ative and positive tests, we conclude that for systems with 
^members > 20, the DS Test can reliably identify real sub- 
structure if the P-values method is employed using a confi- 
dence level of either 0.01 or 0.05. For groups with n mom b cr s < 
20, we find that the percentage of groups with substructure 
identified by the DS Test should be taken as a lower limit. 



3 SUBSTRUCTURE IN THE GEEC GROUPS 
3.1 The GEEC Group Catalog 

The Group Environment and Evolution Collaboration 
(GEEC) group catalog is based on a set of ~200 inter- 
mediate redshift, 0.1 < z < 0.6, galaxy groups initially 
identified in the second Canadian Netwo rk for Observa- 
tional Cosmology (C NOC2) redshift survey jYee et al.ll2000l ; 
ICarlberg et alj200ll ). The CNOC2 survey observed ~4 x 10 4 
galaxies covering four patches, totalling 1.5 deg 2 in area, in 
the U BVRcIc bands down to a limiting magnitude of Rc = 
23.0. Spectra of more than 6000 galaxies were obtained with 
the MOS spectrograph on the Canada-France-Hawaii Tele- 
scope (CFHT), w ith 48 per cent completeness at Rc = 21.5 
l|Yee et alj kOQCh . The GEEC group catalog includes ex- 
tensive follow-up spectroscopy with the Inamori Magellan 
Areal Camera and Spectrograph (IMACS) (Connelly, J. et 
al., submitted) and L ow Dispersion Surv ey Spectrograph 
(LDSS2) on Magellan (|Wilman et alj|2005l ). as well the Fo- 
cal Reducer and low dispersion Spectrograph (FORS2) on 
the Very Large Telescope (VLT) (Connelly, J. et al., sub- 
mitted). We have also obtained multi- wavelength imaging 
data, which includes: X-ray observations with the X-ray 
Multi-Mirror Mission-Ne wton (XMM-Newton) a nd Chan- 
dra X-ray Observatories (|Finoguenov et al] [2009). ultravio- 
let observations wit h Galaxy Evolution Explorer (GALEX) 
l|McGee et al.l20l"lT ). optical observations with the Advanced 
Camer a for Surveys (ACS) on the Hubble Space Telescope 
(HST) (|Wilman et alj|2009l ), infrared observations with the 
Multi-band Imaging Photo meter for Spitzer (MIPS) on 
the Spitzer Space Telescope (|Tvler et al.ll201.lh . and near- 
infrared observations with Isaac Newton Group Red Imag- 
ing Device (INGRI D) on the William Herschel Telescope 
dBalogh et alj 120091), Infrar ed Array Camera (IRAC) on 
Spitzer (|Wilman et al.ll2008h , and the Son of I SAAC (SOFI) 
on t he New Technology Telescope (NTT) (|Balogh et all 
2009). In addition, improved optical imaging was obtained 
in the uqrizBV R I filters from the CFHT Megacam and 
CFH12K imagers (|Balogh et alj|2009h . 

In addition to the follow-up observations of the CNOC2 
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fields, group memb ership has also been redefined by 
IWilman et all (|2005m usi ng more relaxed algori thm param- 
eters than those used bv lCarlberg et aD l|200lf ). The origi- 
nal search parameters were optimized such that the group- 
finding a lgorithm would iden tify dense, virialized groups, 
while the lWilman et al.1 (2005) catalog includes looser group 
populations that cover a wider range of dynamical states. 
For this reason, the GEEC group catalog is ideal for the in- 
vestigation of substructure within groups, since we are not 
restricted to dense group cores and are able to probe the 
surrounding large scale structure. 



3.2 Analysis of the GEEC Groups 

We apply the DS Test, as described in Section 2, to a sub- 
set of the GEEC groups. Although there are roughly 200 
groups in the GEEC catalog, the majority of systems have 
fewer than 10 members. In Section 2.2.3, we determined that 
in order to obtain a reliable measure of substructure, the 
minimum number of member galaxies for the DS Test is 
^members = 20, which leaves us with a subset of 15 groups. 
These groups are amongst the most massive GEEC groups 
with an average velocity dispersion of ~525 km s _ . 

As previously mentioned, due to the relaxed member- 
ship allocation algorithm parameters, some of the GEEC 
groups have relatively large radial extents, and can be larger 
than the suggested max imum group virial radius (r2oo) of 
1.0 Mpc |m 

amo nl l2007D . but for the purposes of detecting 
substructure we elect not to apply any radial cuts to our 
groups for our main analysis. Our definition of substructure 
is liberal, in that we include structure that may be infalling 
or structure that is in nearby large scale structure but not 
necessarily bound to the host group. Thus, to ensure that we 
do not eliminate any possible detection of substructure, we 
analyze all galaxies i dentified as group me mbers by the FOF- 
algorithm applied in lWilman et al.l (|2005l ). Our decision not 
to apply radial cuts is further justified give n that substruc- 
ture i s often in group and cluster ou tskirts l|West fc Bothunl 
ll990l ; IZabludoff fe Mulchaey||l998al ). 

For the distribution parameters we estimate V\ OC!3l and 
z7, from Equation [1] as the group and local (i.e. ith galaxy 
and its ^/n mcm b cr s nearest neighbours) canonical mean ve- 
locity, and o"i ocal and a, also from Equation [T] as the local 
and group intrinsic velo city dispersion, com puted following 
the method outlined in IWilman et al.1 (|2005l l. The disper- 
sion uncerta inties are computed using the jackknife method 
l|Efronll 19821 ) and are given in Tables [2] and [3] 

Of the 15 GEEC groups with n mom b cr s > 20 we find 
that 4 groups, (~27 per cent) are identified as having sub- 
structure at the 99 per cent confidence level (c.L). In Table 
we list group properties, A/n mem bcrs-values, and P- values 
for the groups with substructure, identified as the systems 
that have P- values of less than 1 per cent. In Table [3] we list 
the same values, but for groups without substructure, that 
is systems that have P-values greater than 1 per cent. Also, 
it should be noted that although we list the A/n mcm b C rs 
critical values in Tables and [3] we rely on the P-values 
to identify substructure in our groups. In addition, Tables 
[2] and [3] lists the dynamical properties of the groups to be 
discussed in Sections 4.1 and 4.2. In Fig. [2] we plot both the 
intrinsic velocity dispersion of the group (top) and the total 
number of group members (bottom) versus the mean group 
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Figure 2. Top: Intrinsic velocity dispersion Oi n t versus the group 
redshift for the galaxy groups in our sample. Solid circles indicate 
groups with no substructure and open triangles indicate groups 
with identified substructure according to the DS Test, at a 99 
per cent confidence level. The smaller open circles represent the 
dispersions of all the groups in the GEEC catalogue, the major- 
ity of which are not used in this analysis. Our ri mem b ers > 20 
GEEC sample tend to have higher velocity dispersions. Bottom: 
"members versus the group redshift for the galaxy groups in our 
sample. Symbols are the same as the plot above. 



redshift. Groups with and without substructure both span 
a wide range of velocity dispersions and group membership 
indicating that there is no apparent redshift bias with re- 
gards to group velocity dispersions or n mcm bors for the DS 
Test. 

Although the minimum membership cut needed to de- 
termine a reliable percentage of groups that contain sub- 
structure is n mcm bors = 20, we can still apply the DS Test 
to systems with fewer members and establish a lower limit 
on this fraction. Including groups with as few as ten mem- 
bers increases our sample size to 63 systems. Applying the 
same methodology described above we find that 11 of our 
63 groups (~17 per cent) are identified as having significant 
(> 99 per cent c.L) substructure. 



3.3 GEEC Groups with Substructure 

We examine the GEEC groups classified as having substruc- 
ture by the DS Test in detail to determine if we can iden- 
tify the regions of localized substructure. By examining the 
available position and velocity information, we can find col- 
lections of galaxies that are kinematically distinct from the 
host group. 

In the following analysis, we focus on small subgroups 
of galaxies that may be part of some localized substructure. 
This is done by looking at the Si histogram (Figs. E[a) - 
EJa)), velocity histograms (Figs.[3jb) -E^b)), 'bubble-plots' 
i.e. position plots of the gro up members weighted by exp Si 
Dressier fc Shectmanl 119881 ) (Figs. Etc) - Etc)) and group- 
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Table 2. Group properties and DS statistics for the GEEC groups with substructure 



GEEC Group ID 


^members 


% group 


17 


^/^members 


P value" 


AD Test 


Shape of Velocity 








km s — 1 






Classification 


Dispersion Profile 


25 


28 


0.362 


491±32 


0.693 





non-Gaussian 


rising 


208 


22 


0.269 


530±45 


0.841 





non-Gaussian 


rising 


226 


86 


0.359 


944±17 


1.37 





non-Gaussian 


rising 


320 


29 


0.245 


463±33 


0.751 


0.00383 


Gaussian 


rising 



"Using 100 000 'MC shuffles'. We identify groups with P values < 0.01 as having significant substructure. 



Table 3. Group properties and DS statistics values for the GEEC groups with no substructure 



GEEC Group ID 


"■members 


Z group 




^/^members 


P-value a 


AD Test 


Shape of Velocity 








km s _1 






Classification 


Dispersion Profile 


4 


20 


0.201 


360±39 


0.605 


0.524 


non-Gaussian 


flat 


38 


34 


0.511 


739±21 


0.910 


0.0372 


non-Gaussian 


declining 


104 


27 


0.145 


365±25 


0.659 


0.0501 


Gaussian 


flat 


110 


33 


0.156 


338±22 


0.599 


0.413 


Gaussian 


declining 


117 


27 


0.220 


261±20 


0.412 


0.851 


Gaussian 


flat 


138 


53 


0.438 


743±29 


0.983 


0.354 


non-Gaussian 


flat 


238 


21 


0.408 


606±52 


0.826 


0.135 


non-Gaussian 


flat 


308 


26 


0.224 


511±35 


0.854 


0.610 


Gaussian 


declining 


334 


20 


0.323 


454±30 


0.670 


0.0563 


Gaussian 


flat 


346 


32 


0.373 


439±20 


0.612 


0.232 


non-Gaussian 


flat 


362 


24 


0.460 


666±43 


0.861 


0.0817 


Gaussian 


rising 



"Using 100 000 'MC shuffles'. We only identify groups with P values < 0.01 as having significant substructure. 



centric radial velocity (i.e. C2: mcm b C r — cz SIOup ) weighted po- 
sition plots (Figs. E^d) -Od)). 

The Si and velocity histograms provide an estimate on 
the amount of substructure and the dynamical state of the 
groups. The Si histogram gives an overall view of the kine- 
matic deviations and the velocity histogram can be used 
to identify non-Gaussian features, such as multiple peaks. 
In order to look for local regions of substructure, one must 
look at both the 'bubble- plot' and velocity weighted posi- 
tion plots in Figs. O - [5] The 'bubble-plots' allow for the 
visual identification of candidate regions of local substruc- 
ture. Since the size of the symbols in the 'bubble-plot' scales 
with a galaxy's Si value, large symbols correspond to strong 
local kinematic deviations from the global values. Thus, a 
collection of neighbouring galaxies with similarly large sym- 
bols, such as region A in Fig. [Sfc), could indicate a kine- 
matically distinct system. In order to confirm that these 
candidate regions are truly distinct, one must check that 
the candidate substructure galaxies also have similar veloc- 
ities, since the sign or direction of the galaxy velocity is 
not taken into account in the DS Test (Equation [TJ . To de- 
termine this, we look at the group-centric velocity weighted 
position plot to see if the candidate local substructure galax- 
ies have similar velocities. In Figs. [3jd) to [6jd) the red 
symbols correspond to positive group-centric velocities (i.e. 
cZmember — cZgroup > 0), blue symbols correspond to negative 
group-centric velocities (i.e. cz me mber — cz group < 0) and the 
symbol size scales with the magnitude of the velocity off- 
set. We only identify galaxies as part of local substructure 
if neighbouring galaxies have similar large kinematic devia- 



tions, shown in the 'bubble-plots' and similar group-centric 
radial velocities, shown in the weighted position plots. 

We will now discuss each of the four GEEC Groups with 
substructure in detail. Using the methodology described 
above, we search for candidate local regions of substructure 
in our sample. GEEC Group 25 (Fig. [3]) has a collection of 
five galaxies, just south-east of the group centre, that all 
have high Si values and comparable velocities. Similarly in 
GEEC Group 208 there are seven galaxies that lie south- 
west of the group centre, with equally high Si values (Fig. 
3{c)). Though, when we look at these same galaxies in the 
velocity weighted position plot (Fig.[4fd)), only five of these 
seven galaxies have comparable radial velocities. This result 
highlights the importance of looking at both the 'bubble- 
plot' and velocity weighted positions plots, as galaxies with 
large kinematic deviations may be correlated in position- 
space, but not in velocity-space. 

In Fig. O we show the substructure analysis plots for 
GEEC Group 226 and from the 'bubble- plot' we see that 
there are two possibly distinct regions of localized substruc- 
ture. The first region lies directly north-east of the group 
centre (region A in Fig. [5jc)) and includes the galaxy with 
the highest Si value. All of the members within this given 
substructure have similar radial velocities (Fig. EJd)). The 
second region of interest contains eleven galaxies that lie in 
the very north-east corner of the 'bubble-plot' (region B in 
Fig. EJc)). Again, all of the galaxies within this particular 
substructure have similar group-centric velocities, though in 
the opposite direction of the members in region A. 
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Figure 3. GEEC Group 25. (a): Si histogram, (b): Histogram of 
the velocity distribution, where the dashed-line indic ates the best 
fitting Gaussian velocity dispersion, (c): Dressier & ; Shectmanl 
(1988) 'bubble-plot' where the galaxy symbols scale with exp(<5i). 
(d): Position plot where the galaxy symbols scale with group- 
centric velocity (i.e. ) /350), blue symbols 
correspond to galaxies with negative group-centric velocities (i.e. 

2 g roup < 0) and red symbols correspond to galaxies 



^member C* 

with positive group-centric velocities (i.e. cz mem b er - 



>0). 
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Figure 5. Same as Fig. [3] but for GEEC Group 226 and the ve- 
locity weighted plot now scales as exp (ez mem b e r — ez gr0 up) /600. 
The dotted box encompasses the first identified region of local 
substructure (region A) and the long dashed box encompasses 
the second region of local substructure (region B). 
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Figure 6. Same as Fig. [3] but for GEEC Group 320. 



Figure 4. Same as Fig. [3] but for GEEC Group 208 and the ve- 
locity weighted plot now scales as exp (ez mem b e r — cz group ) /400. 
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From the 'bubble-plot' of GEEC Group 320 (Fig. He)), 
we see two possible regions of substructure; one just north- 
west of centre and another further south-west of centre. 
The structure near the group centre does not appear to 
be very localized, as the velocities, though in the same di- 
rection, have significantly different magnitudes (Fig. EJd)). 
The second region of high 8i values may actually be two 
separate structures. The velocity weighted position plot in 
Fig. EJd) reveals that three of the galaxies in the structure 
have similar negative group-centric velocities (i.e. C0 me mbef — 
cZgroup < 0) and two have similar positive velocities (i.e. 
cZmembor — cz gr0 up > 0). Therefore, we identify the collec- 
tion of three galaxies, with the negative group-centric radial 
velocities, as the best candidate of substructure within this 
system. 

In each of the GEEC groups with substructure, it ap- 
pears that our identified regions of localized substructure lie 
on the outskirts or edges of the group, and may either be 
infalling onto a pre-existing system or in the nearby large 
scale structure, but not necessarily bound or infalling. This 
result is in agre ement with similar studies of substructure 
in local groups (Zabludo ff fc MulchaevMl998aT ) and clusters 
|West fc Bothunlll99d) . ~ 

It should also be noted that although the DS Test has 
been shown to be reliable for group-sized systems, the num- 
ber of members within the substructure can affect the results 
of the statistic. As discussed in Section 2.2, and shown in 
more detail in the Appendix, fewer(more) true members of 
substructure can increase(decrease) the rate of false nega- 
tives. However, we also show that the rate of false positives 
is consistently low for all of our mock groups (see Section 
2.2). Thus, any detection of substructure in these systems is 
likely to be real. 



3.4 Is the localized substructure gravitationally 
bound to the group? 

The local regions of substructure we detect lie on the group 
outskirts and are possibly bound and infalling, or not bound, 
but close in large scale structure. In order to help distinguish 
between these two possibilities, we apply a simple bound 
test to estimate whether or not the regions of localized sub- 
structure are gravitationally bound to the host group. This 
is done by computing the limits for bound sy stems, using 
a vari ation of the virial theorem as discussed in iBeers et al.l 
i| 19821 ) and determ i ning i f the substructure falls within these 
limits. IBeers et alj |l982T ) state that for a two-body system 
on a linear orbit, the Newtonian limit for gravitational bind- 
ing, projected onto the sky is given by 



where, 



VrRp < 2GM sin 2 a cos a, 



V r — Vsina, R p — J? cos a, 



(5) 



(6) 




and where a is the angle between the line joining the 
two-body system and the plane of the sky (see fig. 7 of 
IBeers et al.l 1 198^ ) . M is the total mass of the entire sys- 
tem (substructure plus host group), and R and V are the 
true (3-dimensional) positional and velocity separations be- 



200 400 

V r (km s" 1 ) 

Figure 7. a - V r plot for GEEC Groups (a) 25, (b) 208, (c) 
region A and the host group of 226, (d) region B and the host 
group of 226 (d) and (e) 320 . The unshaded regions correspond 
to bound solutions and shaded regions correspond to unbound 
solutions of the virial theorem given by Equation \E\ The solid 
line corresponds to the measured value of V r (i.e. the line-of- 
sight velocity difference between the substructure and host group 
centres). The dashed lines correspond to one-sigma deviations, 
taken to be the errors on the intrinsic velocity dispersion. 



tween the two objects. V r is the line-of-sight relative velocity 
between the two bodies and R p is the projected separation, 
both of which are measurable quantities. 

The unknown quantity in Equation [5] is the projection 
angle a. Thus, to compute the limit between bound and 
unbound systems, one must work in a — V r space to de- 
termine the probability that the system is gravitationally 
bound for any given projection angle. This is achieved by 
setting 0° < a < 90° and solving for V r in Equation [S] 
producing a distinct line in a — V r space that clearly sepa- 
rates bound and unbound solutions (Fig. 0. One can then 
compute the probability of a bound solution, for a given 
projection angle, using the a — V r plot. 

To apply this methodology to our group sample, we 
treat our identified regions of local substructure as one-body 
and the remaining galaxies as the second-body, which we 
refer to as the 'host' group. The total mass of the system is 
taken to be the virial mass, M200 (i-e. the total mass within a 
radius that encloses a mean density of 200 times the critical 
density of the Universe at the redshift of the galaxy), of 
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the GEEC Groups as computed in iBalogh et all (|2009r ) and 
given by 



Table 4. Groups properties and DS statistics values for the 
GEEC Groups with a 1.0 Mpc radius cut 



M 200 = 



(7) 



G lOHo (1 + z) 15 ' 

where a is the measured intrinsic velocity dispersion. The 
measured V r and R p values are taken to be the distance 
between the P-band luminosity- weighted centres of the local 
substructure and the 'host' group centre, along the line-of- 
sight (V r ) and projected on the sky (-Rp). 

In Fig. [7] we plot the a — V r plots for the five candidate 
regions of local substructure, as discussed in Section 3.3. 
The shaded regions indicate the areas spanned by unbound 
solutions, as given by Equation [5] the solid black vertical 
line is the measured value of V r , and the long-short dashed 
vertical lines indicate one sigma deviations, taken to be the 
intrinsic velocity dispersion error. 

Using the methodology described above and Fig. [7] we 
conclude that the identified local substructures in GEEC 
Groups 25 and 226 are not bound to the host group, while 
for GEEC Groups 208 and 320 the detected substructure is 
likely bound to the host group. 



3.5 Substructure within 1 Mpc of the Group 
Centroid 

In the previous section we analyzed a subset of 15 GEEC 
groups, with n mem bers > 20, without applying any radial 
cuts. Here we apply a 1.0 Mpc cut, which is the s uggested 
maximum virial radius for groups (|Mamon|[2007l ). on the 
same subset of groups and re-apply the DS Test. Applying 
this radial cut, while still requiring a minimum number of 
20 member galaxies, reduces our sample from 15 to 5 groups 
(GEEC Groups 110, 138, 226, 308 and 346). With this radial 
cut, we find that all 5 groups are identified as not having sub- 
structure (see Table S} according to the DS Test. The only 
group that was previously identified as having substructure, 
prior to the 1.0 Mpc cut, is GEEC Group 226. In Fig. [5] it is 
evident that although there are galaxies with relatively high 
5i values close to the group centre, the members with the 
highest Si values lie near the edge of the group. In fact the 
most significant feature in the 'bubble-plot' (Fig. [5| is on 
the top-right corner of the plot, far from the group centre. 

If we include all of the groups with n mcm b crs > 10 after 
a 1.0 Mpc radial cut, we find that 2 out 33 groups (~6 per 
cent) of our sample contains significant substructure. Again, 
we note that for systems with fewer than 20 members, the 
DS Test can only provide a lower limit on the amount of 
substructure present. 



4 CORRELATIONS BETWEEN 
SUBSTRUCTURE AND OTHER 
INDICATORS OF DYNAMICAL STATE 

Having identified the GEEC groups that contain substruc- 
ture, we can also look at other dynamical properties, i.e. 
the velocity distributions and velocity dispersion profiles 
(VDPs), to determine if there are any correlations with sub- 
structure. 



GEEC Group ID 


n a 

members 


V 




P-value 






km s — 1 


km s _1 




110 


26 


-15.3 


350 


0.403 


138 


23 


-226 


730 


0.792 


226 


25 


-174 


847 


0.110 


308 


25 


2.52 


512 


0.507 


346 


26 


-80.4 


434 


0.128 



"Group membership after the 1.0 Mpc radius cut 

4.1 Comparison with the Dynamical 

Classification of Velocity Distribution 

If a correlation does exist between substructure and re- 
cent galaxy accretion, one would expect that the groups 
with substructure should also be dynamically co mplex, with 
perha ps non-Gaussian velocity distributions. In iHou et al.l 
(2009), we established a classification scheme to distinguish 
between dynamically relaxed and complex groups. Using the 
Anderson-Darling (AD) goodness-of-fit test, we are able to 
determine how much a system's velocity distribution de- 
viates from Gaussian. This is done by comparing the cu- 
mulative distribution function (CDF) of the ordered data, 
which in our case is the observed velocity distribution, to 
the model Gaussian empirical distr ibution function (EPF) 
using computing formulas given in iD'Agostino fc Stephens! 
( 1986). The AD statistic is then converted into a probability, 
or P- value, us i ng res ults determined via Monte Carlo meth- 
ods in iNelsonl l| 19981 ). A system is then considered to have a 
non-Gaussian velocity distribution if its computed P-value 
is less then 0.01 corresponding to a 99 per cent c.l. 

We now apply this scheme to our sample of 15 GEEC 
groups, with n mem b er s > 20 and no radial cut, to compare 
the dynamical state with the detection of substructure The 
results of our dynamical classification scheme indicate that 
8 of the 15 groups are classified as having non-Gaussian ve- 
locity distributions, at the 99 per cent c.l., and are thus dy- 
namically complex (see Tables [2] and [3} . Of the four GEEC 
groups with substructure only GEEC Group 320 shows a 
velocity distribution consistent with being Gaussian. Five 
groups, GEEC Groups 4, 38, 138, 238, and 346, are identi- 
fied as being dynamically complex, but do not contain any 
subs tructure accord i ng to the DS Test. 

IPinknev et al.l (|l996t ) found that different statistical 
tools (i.e. 1-D, 2-D and 3-D tests) probe the dynamical 
state of a system at varying epochs. Using iV-body simula- 
tions, these authors determined that 1-D tests, such as the 
AD Test are most sensitive to scenarios when substructure 
passes through the core of the host group (i.e. core-crossing). 
During this time substructure can become spatially mixed 
within the host group, and if the substructure is loosely 
bound then it may be difficult to detect with 3-D tests, such 



3 The dynam i cal cl assifications in this paper differ from those 
in IHou et al. The reason for this difference is that in 

IHou et al.l ||2009| 1 we applied a 1.0 Mpc radius cut to the GEEC 
groups, whereas no radial cut is applied in this analysis. 
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as the DS Test. Thus, groups with non-Gaussian velocity dis- 
tributions may contain substructure that is missed by the 
DS Test. 



4.2 Comparison with the Velocity Dispersion 
Profiles 

In a study of th e VD Ps of galaxy clusters, 
M cnci fc Fusco-Femiand (|l996h found a correlation be- 
tween the efficiency of merger activity in the cluster core 
and the shape of the VDP, where actively interacting 
systems had strongly rising pro files. We presente d a similar 
correlation for galaxy groups in lHou e~ al. (2009), where we 
found that groups classified as dynamically complex (i.e. 
non-Gaussian velocity distributions) also had rising VDPs. 
We now look at the VDPs of the current GEEC group 
sample to see if there is a similar relationship between the 
shape of the profile and the detection of substructure. 

The VDPs are com puted following the method outlined 
in lBergond et al] (|2006T ). Unlike traditional methods of com- 
puting binned projected velocity dispersions this technique 
generates a smoothed VDP. This is done by using a 'moving 
window' prescription, which takes into account the contri- 
bution of every radial velocity measurement at each com- 
puted radius. The values are binned with an exponentially 
weighted moving window given by 



Wi (R) — — exp 



(R-Ri 



2a% 



(8) 



where a a is the width of the window, which can be constant 
or a function of radius R, and the Ri's axe the radial posi- 
tions of the members of the system. The projected velocity 
dispersions are then defined as: 



a p (R) 



Ej Wi(R)(xi- 



(9) 



where the Xi's are the radial velocities and x is the mean 
velocity of the system. 

We compute VDPs for our GEEC group sample and the 
profiles are shown in Fig. fusing a widow width, <tr, equal 
to one-third the maximum group radius. It should be noted 
that the projected velocity dispersions do not include any 
redshift or instrumental error corrections, so the intrinsic 
dispersion values (Tables [2] and [3} are generally lower than 
the projected velocity dispersions. 

Comparing the profiles of the groups with and without 
substructure, we find that all four GEEC groups identified 
as having substructure (Groups 25, 208, 226 and 320) also 
have strongly rising profiles. In contrast, almost all of the 
groups with no detected substructure, with the exception of 
GEEC Group 362, either have flat or generally decreasing 
VDPs, within the intrinsic velocity dispersion error. See Ta- 
bles [2] and for a description of the shape of the VDP for 
each GEEC group. Thus, we do indeed observe a correlation 
between detectable substructure and rising VDPs. As previ- 
ously mentioned, studies of rich galaxy clusters suggest that 
a strongly increasing profile may be a signature of merger ac - 
tivity or galaxy interactions (|Menci fc Fusco-Femianoll9 96). 
but an alternative explanation of rising VDPs i s the presence 
of subclumps with differe nt mean velocities l|Girardi et al.l 
1 19961 : iBarrena et al.ll2007f ). In our case, it is likely that the 
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Figure 8. Velocity dispersion profiles (VDPs) for the GEEC 
Groups with n > 20: panel (a) Group 4, (b) Group 25, (c) Group 
38, (d) Group 104, (e) Group 110, (f) Group 117, (g) Group 138, 
(h) Group 208, (i) Group 226, (j) Group 238, (k) Group 308, (1) 
Group 320, (m) Group 334, (n) Group 346 and (o) Group 362. 
The plots with asterisks in the top-left corner indicate groups 
that have been identified as having significant substructure. The 
intrinsic velocity dispersions, with errors, for each group can be 
found in Tables [2] and [3] 



increasing profile is being caused by the kinematically dis- 
tinct substructure, which has a different mean velocity from 
the host group. 



5 CORRELATIONS BETWEEN 

SUBSTRUCTURE AND GALAXY 
PROPERTIES 

Thus far, we have looked at possible correlations between 
substructure and other indicators of the dynamical state of a 
system. We now compare properties of the member galaxies 
to see if there are any differences between the galaxies in 
groups with and without substructure. 



5.1 Substructure and Colour 

In the following analysis we compare the 0,4 (g — r) colours, 
which have been corrected for g alactic extinction a nd k- 
corrected to a redshift of z = 0.4 fBa logh et al.l [2009) . and 
blue fractions, /{,, for the groups with and without substruc- 
ture. 

In addition to extinction and k-corrections, we also ap- 
ply a completeness correction to address the differing spec- 
troscopic coverage. We apply magnitude weights that de- 
pend on whether the group had follow-up spectroscopy. 
The weights we ap ply are similar to those derived in 
IWilman et all (|2005h . except we do not include any radial 
weights. We compute weights in r-band magnitude bins of 
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Figure 9. Left: 0,4 (g — r) colour histograms for the GEEC groups 
with substructure (solid line), which has be normalized to match 
the number count of the groups with no substructure and for the 
GEEC groups with no detectable substructure (dashed line). Both 
histograms have r-band magnitude based weights which have been 
computed to take into account the differing spectroscopic cover- 
age between the original CNOC2 survey and the follow-up Mag- 
ellan survey. Right: Same as figure on the left, expect a 1.0 Mpc 
radial cut has been applied to all groups in the sample. 



0.25, and up to a limit of r = 22.0, which is the limit of 
the unbiased Magellan spectroscopy. These weights are then 
applied to all of the member galaxies in our sample. 

In Fig. [5]we show the completeness weighted 0A {g — r)- 
histograms for the galaxies, with r < 22.0, in groups with 
substructure (solid line) and for those in groups with no 
detected substructure (dashed line). From this figure, it is 
clear that both histograms are bimodal with a distinct red 
sequence and blue cloud. Although both colour distribu- 
tions have the expected bimodal shape, it is obvious that 
galaxies in the groups with and without substructure come 
from very different parent colour distributions. The groups 
with no identified substructure have a well populated red se- 
quence, while the groups with substructure appear to have a 
much more dominant blue cloud. Results from a two-sample 
Kolmogorov-Smirnov (KS) Test on the unbinned 0A (g - r) 
distributions show that these two samples are distinct at the 
99 per cent c.l. 

We also plot colour distributions of the groups with 
and without substructure with a 1.0 Mpc radial cut applied 
to all groups in our sample (Fig. [9] right). Although, the 
blue cloud for the groups with substructure in less popu- 
lated when a radial cut is applied, it is still more populated 
than the groups without substructure and the two distribu- 
tions are still statistically distinct. This result indicates that 
the increase in blue galaxies in the groups with substruc- 
ture is not only coming from galaxies at radi i greater than 
1.0 M pc. Our findings are similar to those of iRibeiro et al.l 
i|2010h who showed that there are many more red galaxies 
in dynamically evolved group systems, even out to 4 virial 
radii. 

To obtain a more quantitative comparison we compute 
the fraction of blue galaxies within each sample. We take 
a multi-stage approach to determine the appropriate colour 
cut needed to distinguish between the red sequence and blue 
cloud. First, we apply an initial colour cut of 0A (g — r) = 1.2, 
based on the minimum value in the colour distribution of 
all the member galaxies in our sample, shown in Fig. 1101 
We then determine a linear fit to all of the galaxies with 
0A (g — r) > 1.2 and set the colour cut to be one standard 




Figure 10. Left: Weighted 0A (g — r) histogram of all the member 
galaxies in our group sample. Right: 0A (g — r) versus r colour- 
magnitude diagram of all the galaxies in our sample, with no 
completeness correction. The solid line indicates the colour cut 
used to distinguish between the red sequence and the blue cloud. 



deviation below the fit to the red sequence. The final colour 
cut is determined to be 



0A (g - r) = -0.0236(°- 4 r) + 1.810, 



(10) 



and is represented by the solid line in the colour-magnitude 
diagram shown in Fig. 1101 

The blue fraction, ft, is then computed as the ratio of 
galaxies with 0A {g — r) values that fall below Equation \W\ to 
the total number of galaxies in the sample. The error in the 
blue fraction is computed using con fidence interval s (CIs) 
derived from the beta distribution (ICamerodlioiol ). This 
method has been shown to more accurately determine CIs, 
especially for small samples, over traditional methods such 
as Poisson errors, which systematically under-estimates the 
width of the CIs. We find that f b — 41±| per cent for the 
groups without substructure and ft = 69±| per cent for the 
groups with substructure. Thus, the groups with substruc- 
ture have a significantly higher blue fraction, which is clear 
in the 0,4 (g — r)-histograms of Fig. [9] We note that although 
these blue fractions are derived from weighted colour distri- 
butions, they are in agreement with the unweighted valuesQ 
suggesting that the applied magnitude weights do not affect 
the observed differences in colour. 

In addition, we also compare the blue fractions of the 
individual galaxies identified as being part of the local sub- 
structure (see Section 3.3) to the other members of the 
group. We find that for the galaxies in the identified regions 
of localized substructure, ft — 74±f per cent, and for the 
galaxies not in the substructure, ft ~ 50±| per cent. This 
suggests that the observed increase in the blue fraction of 
groups with substructure is being enhanced by the galaxies 
in the identified regions of local substructure. 

We now compare our groups with substructure sample 
to a sample of intermediate redshift field galaxies. In Fig. 1111 
we reprodu ce the A (g — r)-his togram for the field galaxies 
in Fig. 6 of B alogh et al.l (|2009l ) (dashed line) , summing up 
all of the counts in each of the quoted Mk b bins in order 
to get the total colour distribution^. We also over-plot the 



= 44 ± 4 per cent for the groups without substructure and 
ft = 68±g per cent for the groups with substructure 
5 It should be noted that although only a fraction of our sample 
actually have measured Mk s values, tho se that do span the e ntire 
range of magnitudes quoted in Fig. 6 of I Balogh et al.l (120091 ') 
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Figure 11. °' 4 (g — r) colour histograms for the GEEC groups 
with substructure (solid line) and for all of field galax ies in the 
GEEC groups found in Fig. 6 of iBalogh et al.l feOOSl) (dashed 
line). It should be noted that the colour distribution for the field 
g alaxies is a sum of all the counts in each magnitude bin of Fig. 6 
of lBalogh et ah 1 1120091 1. since our group galaxy sample covers the 
entire magnitude range. In addition, the group sample has been 
normalized to match the number count of the field sample. 



colour histogram for the galaxies in groups with substruc- 
ture, except we now apply an Mj magnitude cut in order 
to match the Mk b ran ge used in IBalogh et all l|2009l ). Both 
the field galaxies in the Balogh et al. (2009) sample and the 
galaxies in our groups with substructure lack a prominent 
red sequence, and have well populated blue clouds. Despite 
subtle differences in the two colour histograms, a two-sample 
KS Test shows that the 0,4 (g — r) colours of the field galaxies 
and groups with substructure galaxies very likely come from 
the same parent distribution (P- value = 0.67). 

Our blue fractions can be co mpared to those comp uted 
in the zCOSMOS survey, where llovino fc et all (|2010l ) de- 
termined the blue fractions of isolated and group galaxies at 
various redshif t s in t heir sample. At a redshift of z = 0.4, 
llovino fc et all (|20ld ) found f b ~ 70 per cent for isolated 
galaxies and ft ~ 45 per cent for group galaxie^f). From 
these results it is clear that our observed blue fraction of 69 
per cent for the galaxies in groups with substructure is signif- 
icantly higher than the observed zCOSMOS group sample, 
but is in agreement with their field sample. On other hand, 
the fb values for our GEEC groups with no substructure 
roughly agree with the zCOSMOS group sample. 



5.2 Substructure as a Function of Colour, Stellar 
Mass and Star Formation Rates 

In the previous sections we compared the 0A (g — r) colours of 
the galaxies in groups with substructure to those in groups 



6 We look at the Sample II of llovino fc et al l i201Cft . as this 
the data-set that corresponds best to our GEEC sample. 



with no detected substructure. Here we again compare the 
colours of the member galaxies but now as a function of stel- 
lar mass and specific star formation rate (defined as the ratio 
of the star formation rate to the stellar mass - SSFR). Addi- 
tionally, we compare the SSFR distributions of the galaxies 
in groups with and without substructure. It should be noted 
that since only 3 of the 4 original CNOC2 fields were tar- 
geted with GALEX, the sample used in the following anal- 
ysis contains fewer galaxies than used in the colour anal- 
ysis of Section 4.3. In our sample of 15 (n mcm bors > 20) 
GEEC groups, 275 group galaxies have measured SFRs, stel- 
lar masses and colours, while 401 group members have well 
determined colours. 

The stellar masses and star formation rates (SFRs) 
for the GEEC sample were obtained from spectral energy 
distribution (SED) template-fitting to all of the available 
photomet ry. Detailed di s cussio n of the methodology is pre- 
sented in iMcGee et all (|201ll ). but we give a brief sum- 
mary here. The photometry used in the SED-fitting process 
typically included ph otometry in the K , i, r, g, u, NUV 
and FUV bands (see IBalogh et al.ll2009l. for details of the 
GEEC photometry). The observed photometry was then 
co mpared to a large gr i d of m odel SEDs constructed using 
the iBruzual fc Charlotl i|2( )().'> i stellar population synthesis 
code and assuming a C habrier initial mass fu n ction. Follow- 
ing th e methodology of ISalim fc et all i|2007l) . IMcGee et all 
(|2011T) created a grid of models that uniformly sampled 
the allowed parameters of formation time, galaxy metallic- 
ity, an d the two-component dust model of ICharlot fc Fall 
(2000). The star formation history was modelled as an expo- 
nentially declining base rate with random bursts of star for- 
mation of varying duration and relative strength. The model 
magnitudes were obtained by convolving these model SEDs 
with the observed photometric bandpasses at nine redshifts 
between 0.25 and 0.60. ^-minimization was then performed 
by summing over all of the models and taking into account 
the observed uncertainty on each point. The one sigma un- 
certainties in stellar mass, when compared to both mock 
groups and other independent estimates, are on the order of 
0.15 dex and the SFRs h ave been averaged over the last 100 
Myr l|McGee et al.ll201ll ). It should be noted that there may 
be additional systematic uncertainties due to, for example, 
the Initial Mass Function assumed in the fitting procedure. 

In Fig. [12] we show 0A (g - r) versus stellar mass for the 
field galaxies (black dots) , galaxies in groups with substruc- 
ture (closed blue circles), galaxies in groups with no identi- 
fied substructure (closed magenta circles) and the galaxies 
identified as being part of local substructure (open green 
squares). The dotted red vertical line in Fig. 1121 represents 
the stellar mass limit of 1.4 xlO lo M0 at the median red- 
shift (z ~ 0.3) in the GEEC group sample (|McGee et all 
120111 ). It should be noted that this stellar mass limit shifts to 
lower(higher) masses for groups at lower(higher) redshifts. 
From Fig. 1121 it is clear that the galaxies in groups with 
substructure lie preferentially along the blue cloud, as dis- 
cussed in Section 4.3, but we also see that these galaxies 
span a similar mass range as the galaxies in groups with no 
identified substructure. A two-sample KS Test of the stellar 
mass distributions of the galaxies with and without sub- 
structure shows that the two distributions likely come from 
the same parent distribution (P-value = 0.55) The similar 
stellar mass distribution of the two populations tells us that 
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that populate each region of the colour-SSFR space in Fig. 
[T3 with errors com puted using the methodology described 



Cameron] (|2010l ) . Similarly, Table [5] lists the same infor- 
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Figure 12. °- 4 (g — r) versus stellar mass for field galaxies (black 
dots), galaxies in groups with substructure (blue circles), galaxies 
in groups with no substructure (magenta circles) and galaxies 
identified as part of localized substructure (green squares). The 
red dotted line indicates the stellar mass limit (M stc iiar = 1.4 X 
1O 1O M ) at median redshift (z ~ 0.3) of the GEEC sample. 



the DS Test does not simply detect substructure in groups 
which are probed further down the stellar mass function. 

Now we examine the S SFRs for grou p s with and with- 
out substructure. Following iMcGee et alj l|201lf) . we define 
actively star- forming galaxies to have log 10 (SFFR) > — 11 
and the fraction of actively star-forming galaxies to be 

/ 'active — ^-actively star-forming galaxies /^te-tal ■ ^Ve find that for 

galaxies in groups with substructure /active = 63 ±8 per cent 
and for groups with no detected substructure /active = 49±5 
per cent. These active fractions agree with the blue fractions 
found in Section 5.1. 

Both blue and active fractions are used as independent 
indicators of quiescent versus actively star-forming galax- 
ies. However, colour and SSFR probe significantly differ- 
ent timescales, and might be telling us something different 
about the star formation history of galaxies. For instance, 
a dust enshrouded star-forming galaxy would be classified 
as 'red', and therefore quiescent, based on colour but would 
be classified as actively star-forming based on SSFR. A bet- 
ter approach is to look at colour and SSFR simultaneously 
l|Weinmann et al.ll2006t ). We follow this approach in Fig. [T3l 
where we plot 0,4 (g — r) versus SSFR for all the galaxies 
in our sample (top) and with a M st ciiar > 1.4 x 1O 1O M0 
cut applied to the sample (bottom). The colour scheme 
is the same as in Fig. 1121 except the red dotted line now 
corresponds to the division between actively star-forming 
and quiescent/passive galaxies. From Fig. 1131 we can see 
that there is a correlation between colour and SSFR with 
two well populated regions of the plot that correspond to 
'red and passive' and 'blue and active' galaxies, where ac- 
tive refers strictly to actively star-forming galaxies. In Ta- 
ble [5] we list the percentage of all galaxies in our sample 



mation but for galaxies above the stellar mass complete- 
ness limit of 1.4 x 1O 1O M0. From these tables we see that 
the galaxies in groups with substructure have significantly 
more blue and actively star-forming galaxies than groups 
with no substructure, though slightly less than the fraction 
observed in the field. This result indicates that environmen- 
tally driven mechanisms of star-formation quenching are not 
as efficient in groups with observed substructure. We note 
that although the percentages within each region of colour- 
SSFR space differ between the whole versus stellar mass lim- 
ited sample, the general trends remain the same. Since the 
mass-selected sample only includes galaxies with the highest 
stellar masses, we expect a decrease in the 'blue and active' 
region due to stellar mass trends (i.e. high er mass galax- 
ies are preferentially mor e red and passive) l|lovino fc et al.l 
|2010| ; iPeng fc et al.llioToh . 



While the majority of galaxies are either 'red and pas- 
sive' or 'blue and active', a non-negligible fraction appear 
to lie in the other two regions of Fig. 1131 Even if we take 
into account the small uncertainties in c olour (typically 0.02 
mags) and the uncertainties in SSFR l|McGee et alj|201ll . 
quote one sigma errors on the order of 0.25 dex), a measur- 
able fraction of galaxies still remain in these two regions. 
From Table [5] we see that ~9 per cent of all galaxies in the 
field and group samples reside in the 'blue and passive' re- 
gion. This value is hi gher than the ~1.1 per cent observed 
by Wcinmann et all il2006l) , thou gh appears similar to the 



results of lLara-Lopez et alj (|2010l ). Note that both these re- 
sults are based on SDSS galaxies. Although we do observe 
a substantial population of 'blue and passive' galaxies in 
Fig. 1131 we acknowledge that obtaining accurate measures of 
SSFR for ga laxies with low SFR s is notoriously difficult. As 
discussed in lMcGee et al.l (|201ll ) the errors in the measured 
SSFRs increase for galaxies with lower values and some of 
these galaxies may have underestimated SSFRs. 



The final region of the 0A (g - r)-SSFR plot in Fig. [13] 
corresponds to 'red and active' galaxies. While this popu- 
lation is negligible for the groups with substructure sam- 
ple and small for the field sample (see Table [5}, it contains 
~11±3 per cent of the galaxies in the groups with no de- 
tected substructure and ~14±| for the stellar mass limited 
sample. Although this population may be a result of dusty 
star-forming galaxies or edge-on discs with strong extinc- 
tion, it seems unlikely that these galaxies would preferen- 
tially be found in groups with no substructure. An alterna- 
tive explanation is that thi s population could be the 'tran- 
sition' galaxies observed by I Wold (|200gT> . These galaxies are 
still star-forming but at a lower rate (~4 times) than the field 
and have more obscured star formation resulting in weak op- 
tical signatures (i.e. not blue). A population of 'transition' 
galaxies could explain why there are more 'red and active' 
galaxies in groups with no substructure, as a relaxed system 
could contain galaxies that are being quenched but have not 
had their star formation completely cut off. 
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Table 5. Percentage of all galaxies in our sample within a given region of 0,4 (g — r) versus SSFR space 



°' 4 (g-r)-SSFR region 


Galaxies in groups with 
substructure 


Galaxies in groups with 
no substructure 


Field galaxies 


red and passive" 


28 ±6 


42±1 


17 ±1 


red and active 6 


3±| 


11±£ 


4± 1 


blue and passive 


9±3 


9±f 


9±? 


blue and active 


60±| 


38±§ 


70±2 



a Passive denotes quiescent galaxies with log 10 (SFFR) < —11 yr 
b Active refers to actively star- forming galaxies with log 10 (SFFR) > — 11 yr — 



Table 6. Same as Table \E\ except for galaxies above the stellar mass completeness limit of M s t c ii ar > 1.4 X 1O 1O M0 



0,4 (9 - O-SSFR region 


Galaxies in groups with 


Galaxies in groups with 


Field galaxies 




substructure 


no substructure 




red and passive 


42 ± 8 


57 ±7 


33 ± 2 


red and active 


4±| 


14±| 


8±1 


blue and passive 


13±1 


8 ±3 


11 ± 1 


blue and active 


41±8 


21±g 


48 ± 2 



5.3 Implications of the Observed Properties of 
Groups with Substructure 

The field-like colour distribution of the groups with sub- 
structure and the fact that the substructure galaxies are 
found in the group outskirts may have implications for the 
nature of environmental effects in galaxy evolution. It is well 
known that the properties of galaxies depend, at least in 
some part, on their local environment (|Postman fc Gellerl 
Il984l : iDressler et all Il997f ). Galaxies that reside in dense 
environments, such as groups or clusters, generally ex- 
perience some form of star-formation atten uation due to 



19721; lAbadi et all 1 19991; iQuilis et al. 


2000), strangulation 


dLarson et al.lll98d; Baloeh et al.ll200C 


Kawata & Mulchaevl 


20081) or mergers and interactions (|Toomre & Toomrell 19721; 


Broufth et al.l 120061). However, the precise details of the 



galaxy transformation process (i.e. exactly when and where 
quenching occurs, which mechanisms dominate in the differ- 
ent environment, etc.) are still unclear. 

Our substructure analysis in the GEEC groups sug- 
gests that the identified local substructure in our sample, 
which in some cases appears to be infalling, do not feel any 
strong environmental effects from the host group. The ob- 
served colours and SSFR's of the galaxies in groups with 
substructure are significantly more blue, active and remark- 
ably field-like when compared to galaxies in groups with 
no substructure. This suggests that an enhanced fraction of 
red galaxies - either/both via the suppression of star for- 
mation or/and dust obscuration - only happens in relaxed 
groups with no detected substructure. Thus, any environ- 
mental effects felt by infalling substructure galaxies do not 
likely occur until well inside the group potential. 

Recent studies of star-formation as a function of 
group- or cluster-centric radius have produced conflict- 
ing results with regards to the radius at which environ- 
mental effects become observable. Similar to our results, 



IWetzel et al.l (|201ll ) conclude that galaxies do not show sup- 
pressed star formation outsi de the virial radius. In contrast, 
Ivon der Linden et alj l|2010l ) state that suppressed star for- 
mation could be detected in SDSS clusters out to ~ several 
virial radii. Such differences are likely sensitive to group or 
cluster finding algorithms as well as membership assignment. 

In a study of the effects of environment on the colours 
of galaxies in the SDSS survey, I Wilman et all (|2010h take a 
different approach to classifying environment. Rather than 
using a catalog derived from a group- finding algorithm, these 
authors parametrize the environment using non-overlapping 
annular measurements of density on independent scales, al- 
lowing for comparison of envir onmental effect s at v arious 
radii. Based on their analysis, I Wilman et all (|2010l ) con- 
cluded that the fraction of red galaxies correlated with local 
density only up to scales o f ~1 Mpc, which is similar to our 
results, as well as those of IWetzel et alj (|201ll ). Though nu- 
merous and independent analyses of environmental effects 
on galaxy evolution seem to indicate that star-formation 
truncation does not occur until galaxies are well-inside the 
group/cluster potential, there are studies that suggest the 
contrary. Clearly more work, both observational and the- 
oretical, is needed to better understand when and where 
environmental effects become observable. 



6 CONCLUSIONS 

We have studied the Dressler-Shectman Test for substruc- 
ture to determine the sensitivity and reliability of this test 
for group-sized systems. Using mock groups with and with- 
out substructure, generated using Monte Carlo methods, we 
find that the DS Test can reliably be applied to groups 
with more than 20 members, if the probabilities, or P-value, 
method is used with a high confidence level of 99 or 95 per 
cent. We also find that for groups with 10 < n mcmbc rs < 20, 
the DS Test cannot detect all of the substructure within a 
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Figure 13. Top: 0A (g - r) versus specific star formation rate 
(SSFR) for field galaxies (black dots), galaxies in groups with 
substructure (blue circles), galaxies in groups with no substruc- 
ture (magenta circles) and galaxies identified as part of localized 
substructure (green squares). The red dotted line at log 10 (SFFR) 
yr _1 indicates the division between active (log 10 (SFFR) > —11 
yr _1 ) and passi ve (Iog ln (SF FR) < —11 yr -1 ) galaxies deter- 
mined in lMcGee et al.l fooilh . Bottom: Same as figure on the left 
except with a M ste n ar > 1.4 X 1O 1O M0 stellar mass cut applied. 



system, but it can be used to determine a reliable lower limit 
on the amount of substructure. 

Of the 15 rich GEEC groups, with a velocity disper- 
sion range of ~260-950 km s _1 and 20 < n me mbers < 90, 
we find that 4 groups are identified as having significant 
substructure. Further analysis indicates that 2 of these sys- 
tems, GEEC Groups 208 and 320, likely have gravitationally 



bound local substructure that lies on the group outskirts and 
could be accreting onto the system. 

We then looked at various dynamical and galaxy prop- 
erties to search for correlations with the presence of sub- 
structure. The main results of this analysis are; 

(i) The majority of groups with detected substructure 
also have non-Gaussian velocity distributions; 

(ii) The shape of a group's velocity dispersion profile 
(VDP) correlates with the detection of substructure, where 
GEEC groups with substructure have rising profiles; 



(iii) The' 



colour distributions of the groups with 



and without substructure are found to be significantly dif- 
ferent, and the colour distribution for the galaxies in groups 
with substructure is similar to the field distribution; 

(iv) Groups with substructure have a significantly higher 
fraction of blue galaxies, as do the galaxies within identified 
regions of localized substructure; 

(v) Groups with substructure have a larger fraction of 
actively star-forming galaxies (log 10 (SFFR) > —11 yr _1 ), 
when compared to groups with no identified substructure; 

(vi) There is a measurable fraction of galaxies that popu- 
late the 'red and active' region of 0,4 (<? — r)-SSFR space and 
we find that this fraction is significantly higher in groups 
with no substructure for both the whole and stellar mass 
limited samples. 

In conclusion, we find that a considerable fraction of 
intermediate redshift galaxy groups contain significant sub- 
structure, which suggests that like massive clusters, groups 
grow hierarchically through the accretion of smaller struc- 
tures. The field-like colour distribution and measured SSFRs 
of the galaxies in groups with substructure, combined with 
the location of the substructure, suggests that these galax- 
ies are not experiencing any form of environment-related 
star-formation quenching. To fully understand the results 
presented here within the context of galaxy evolution will 
require the use of sophisticated modelling. To this end we 
plan to duplicate this analysis on a sample of semi-analytic 
group s obtained from GALFORM simulations (I Bower et al.l 
2006) and compare these with our observational results. 
With this we hope to be able to better understand the na- 
ture of substructure identified by the DS Test. 



7 ACKNOWLEDGMENTS 

We would like to thank the CNOC2 team for the use of 
unpublished redshifts. A.H, L.C.P, and W.E.H would like to 
thank the National Science and Engineering Research Coun- 
cil of Canada (NSERC) for funding. 



REFERENCES 

G., Moore B., Bower R. G 



1999, MNRAS, 308, 



Abadi M. 
947 

Balogh M. L., McGee S. L., Wilman D., Bower R. G., Hau 
G., Morris S. L., Mulchaey J. S., Oemler Jr. A., Parker 
L., Gwyn S., 2009, MNRAS, 398, 754 

Balogh M. L., McGee S. L., Wilman D. J., Finoguenov A., 
Parker L. C, Connelly J. L., Mulchaey J. S., Bower R. G., 
Tanaka M., Giodini S., 2011, MNRAS, 412, 2303 



Substructure in the Most Massive GEEC Groups: Field-like Populations in Dynamically Active Groups 



Balogh M. L., Navarro J. F., Morris S. L., 2000, ApJ, 540, 
113 

Barrena R., Boschin W., Girardi M., Spolaor M., 2007, 

A&A, 469, 861 
Beers T. C, Geller M. J., Huchra J. P., 1982, ApJ, 257, 23 
Bergond G., Zepf S. E., Romanowsky A. J., Sharpies R. M., 

Rhode K. L., 2006, A&A, 448, 155 
Berlind A. A., et al. 2006, ApJS, 167, 1 
Bird C, 1994a, ApJ, 422, 480 
Bird C. M., 1994b, AJ, 107, 1637 

Bohringer H., Pratt G. W., Arnaud M., Borgani S., Croston 

J. H., Ponman T. J., Ameglio S., Temple R. F., Dolag K., 

2010, A&A, 514, A32+ 
Bower R. G., Benson A. J., Malbon R., Helly J. C, Frenk 

C. S, Baugh C. M., Cole S., Lacey C. G., 2006, MNRAS, 

370, 645 

Brough S., Forbes D. A., Kilborn V. A., Couch W., 2006, 

MNRAS, 370, 1223 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Burgett W. S., et al. 2004, MNRAS, 352, 605 
Cameron E., 2010, ArXiv e-prints 

Carlberg R. G., Yee H. K. C, Morris S. L., Lin H., Hall 

P. B., Patton D. R., Sawicki M., Shepherd C. W., 2001, 

ApJ, 552, 427 
Chariot S., Fall S. M., 2000, ApJ, 539, 718 
Colless M., Dunn A. M., 1996, ApJ, 458, 435 
D'Agostino R., Stephens M., 1986, Goodness-of-fit Tech- 
niques. Marcel Dekker Inc. 
Diemand J., Kuhlen M., Madau P., Zemp M., Moore B., 

Potter D., Stadel J., 2008, Nature, 454, 735 
Dressier A., Oemler Jr. A., Couch W. J., Smail I., Ellis 

R. S., Barger A., Butcher H., Poggianti B. M., Sharpies 

R. M., 1997, ApJ, 490, 577 
Dressier A., Shectman S. A., 1988, AJ, 95, 985 
Efron B., 1982, The Jackknife, the Bootstrap and other 

resampling plans 
Eke V. R., Baugh C. M., Cole S., Frenk C. S., King H. M., 

Peacock J. A., 2005, MNRAS, 362, 1233 
Finoguenov A., Connelly J. L., Parker L. C, Wilman D. J., 

Mulchaey J. S., Saglia R. P., Balogh M. L., Bower R. G., 

McGee S. L., 2009, ApJ, 704, 564 
Firth P., Evstigneeva E. A., Jones J. B., Drinkwater M. J., 

Phillipps S., Gregg M. D., 2006, MNRAS, 372, 1856 
Geller M. J., Huchra J. P., 1983, ApJS, 52, 61 
Girardi M., Demarco R., Rosati P., Borgani S., 2005, A&A, 

442, 29 

Girardi M., Fadda D., Giuricin G., Mardirossian F., 
Mezzetti M., Biviano A., 1996, ApJ, 457, 61 

Gunn J. E., Gott III J. R., 1972, ApJ, 176, 1 

Hou A., Parker L. C, Harris W. E., Wilman D. J., 2009, 
ApJ, 702, 1199 

Iovino A., et al. 2010, A&A, 509, A40+ 

Katz N., White S. D. M., 1993, ApJ, 412, 455 

Kawata D., Mulchaey J. S., 2008, ApJ, 672, L103 

Knebe A., Miiller V., 2000, A&A, 354, 761 

Lacey C, Cole S., 1993, MNRAS, 262, 627 

Lara-Lopez M. A., Bongiovanni A., Cepa J., Perez Garcia 
A. M., Sanchez-Portal M., Castaneda H. O., Fernandez 
Lorenzo M., Povic M., 2010, A&A, 519, A31+ 

Larson R. B., Tinsley B. M., Caldwell C. N., 1980, ApJ, 
237, 692 

Mamon G. A., 2007, in Saviane I., Ivanov V. D., Borissova 



J., eds, Groups of Galaxies in the Nearby Universe The 
Evolution of Galaxy Groups and of Galaxies Therein, pp 
203-+ 

McGee S. L., Balogh M. L., Wilman D. J., Bower R. G., 
Mulchaey J. S., Parker L. C, Oemler A., 2011, MNRAS, 
413, 996 

Menci N., Fusco-Femiano R., 1996, ApJ, 472, 46 
Moore B., Katz N., Lake G., 1996, ApJ, 457, 455 
Nelson L., 1998, Journal of Quality Technology, 30, 298 
Peng Y.-j., et al. 2010, ApJ, 721, 193 

Pinkney J., Roettiger K., Burns J. O., Bird C. M., 1996, 

ApJS, 104, 1 
Postman M., Geller M. J., 1984, ApJ, 281, 95 
Press W. H., Schechter P., 1974, ApJ, 187, 425 
Quilis V., Moore B., Bower R., 2000, Science, 288, 1617 
Ramella M., Biviano A., Pisani A., Varela J., Bettoni 

D. , Couch W. J., D'Onofrio M., Dressier A., Fasano G., 
Kj0rgaard P., Moles M., Pignatelli E., Poggianti B. M., 

2007, A&A, 470, 39 

Ribeiro A. L. B., Lopes P. A. A., Trevisan M., 2010, MN- 
RAS, 409, L124 
Salim S., et al. 2007, ApJS, 173, 267 

Silverman B. W., 1986, Density estimation for statistics 

and data analysis 
Solanes J. M., Salvador-Sole E., Gonzalez-Casado G., 1999, 

A&A, 343, 733 
Springel V., Wang J., Vogelsberger M., Ludlow A., Jenkins 

A., Helmi A., Navarro J. F., Frenk C. S., White S. D. M., 

2008, MNRAS, 391, 1685 

Springel V., White S. D. M., Jenkins A., Frenk C. S., 
Yoshida N., Gao L., Navarro J., Thacker R., Croton D., 
Helly J., Peacock J. A., Cole S., Thomas P., Couchman 
H., Evrard A., Colberg J., Pearce F., 2005, Nature, 435, 
629 

Summers F. J., Davis M., Evrard A. E., 1995, ApJ, 454, 1 
Taylor J. E., Babul A., 2004, MNRAS, 348, 811 
Toomre A., Toomre J., 1972, ApJ, 178, 623 
Tyler K. D., Rieke G. H., Wilman D. J., McGee S. L., 

Bower R. G., Bai L., Mulchaey J. S., Parker L. C, Shi Y., 

Pierini D., 2011, ApJ, 738, 56 
von der Linden A., Wild V., Kauffmann G., White S. D. M., 

Weinmann S., 2010, MNRAS, 404, 1231 
Weinmann S. M., van den Bosch F. C, Yang X., Mo H. J., 

2006, MNRAS, 366, 2 
West M. J., Bothun G. D., 1990, ApJ, 350, 36 
Wetzel A. R., Tinker J. L., Conroy C, 2011, ArXiv e-prints 
Wilman D. J., Balogh M. L., Bower R. G., Mulchaey J. S., 

Oemler A., Carlberg R. G., Morris S. L., Whitaker R. J., 

2005, MNRAS, 358, 71 
Wilman D. J., Oemler A., Mulchaey J. S., McGee S. L., 

Balogh M. L., Bower R. G., 2009, ApJ, 692, 298 
Wilman D. J., Pierini D., Tyler K., McGee S. L., Oemler 

Jr. A., Morris S. L., Balogh M. L., Bower R. G., Mulchaey 

J. S., 2008, ApJ, 680, 1009 
Wilman D. J., Zibetti S., Budavari T., 2010, MNRAS, 406, 

1701 

Wolf C. e. a., 2009, MNRAS, 393, 1302 

Yang X., Mo H. J., van den Bosch F. C, Pasquali A., Li 

C, Barden M., 2007, ApJ, 671, 153 
Yee H. K. C, Morris S. L., Lin H., Carlberg R. G., Hall 

P. B., Sawicki M., Patton D. R., Wirth G. D., Ellingson 

E. , Shepherd C. W., 2000, ApJS, 129, 475 



18 Hou et al. 

Zabludoff A. I., Mulchaey J. S., 1998a, ApJ, 498, L5+ 
Zabludoff A. I., Mulchaey J. S., 1998b, ApJ, 496, 39 



Substructure in the Most Massive GEEC Groups: Field-like Populations in Dynamically Active Groups 



APPENDIX A: FALSE NEGATIVE RATES IN 
MORE DETAIL 

In Section 2.2.3, we presented the main results of the effects 
of changing various input parameters in our mock groups. 
Here we present tables detailing the specific false negative 
rates obtained and also discuss each free parameter in detail. 
As discussed in Section 2.2.3, we determine the sensitivity 
of the DS Test to each of the free parameters in our mock 

groups (<Tp osition j ^position j ^rcdshift j ^rcdshift 5 ^sub 

and (Jhost) 

by beginning with a 'base' mock group (Table[TJ), which has 
a false negative rate of zero per cent. We then change only 
one parameter at a time to ensure that any change in the 
false negative rate can be directly attributed to the altered 
free parameter. 



Al Angular Size of the Substructure (imposition) 

In Table [ATI we present the false negative rates (listed as a 
percentage) for mock groups with n mern t er values of 10, 15, 
20 and 50 as a function of projected angular size on the sky 
(o"poaition). The first line in the table indicate the results for 
our 'base' groups, which have P-values of either or 1 per 
cent. We then increase the value of ^position and from Table 
I Al I we can see that the DS Test reliably identifies substruc- 
ture with a projected dispersion of up to 0.1 Mpc for groups 
with 10 members, and as large as 0.1 Mpc for groups with 
more then 15 members. If we increase exposition to 0.2 Mpc, 
we find that for groups with roughly 20 members or less, the 
false negative rates increases dramatically, but still remain 
very low (1 per cent) for richer groups with 50 members. 
The general conclusion from Table |A"T1 is that the DS Test 
is not very sensitive to the size of the substructure and that 
even for small groups it can identify real substructure that 
is relatively large (~ 0.1 Mpc). 



A2 Location of Substructure in Position Space 

(fi position ) 

In Table lA"2l we list the false negative rates for mock groups 
with n mem ter values of 10, 15, 20 and 50 as a function of the 
projected radial distance of the substructure with respect to 
the group centroid (e pos ition). It is clear from this Table that 
the DS Test is quite insensitive to e p0 sition- In other words, 
the test can reliably identify substructure that is 'close' to 
the projected group centre, and easily detects structure that 
lies on the group outskirts. 



A3 Velocity Dispersion of the Input Substructure 

(Oredshift) 

In Table lA"3l we list the false negative rates for mock groups 
with n mem ter values of 10, 15, 20 and 50 as a function of the 
velocity dispersion of the input substructure (uredshift)- From 
Table lA"3l it is evident that the DS Test is also insensitive to 
this parameter and can reliably detect substructure with a 
wide range of velocity dispersions for all values of n m cmbcrs- 
Only for groups with a fewer than 20 members and a very 
large dispersion value of 450 km s~ do the false negative 
rates go above 5 per cent. 



A4 Location of Substructure Along the 
Line-of-Sight (e reds hift) 

In Table lA"3l we list the false negative rates for mock groups 
with n me mber values of 10, 15, 20 and 50 as a function of the 
location of the substructure along the line-of-sight (e re dshift)- 
This parameter is taken to be a displacement (in km s^ 1 ) 
of the peak in the velocity distribution of the substructure 
with respect to the peak of the host's distribution. Unlike the 
previously discussed parameters, the DS Test appears to be 
extremely sensitive to the location of the substructure along 
the line-of-sight, or more specifically the separation between 
the main groups velocity distribution and the substructure 
velocity distributionty. 

The first set of entries in Table IA3I are the false neg- 
ative rates from mock groups with cr n ost values set to the 
average dispersion of observed GEEC groups with similar 
group membership (Table [TJ. It is clear that if the peaks of 
the host and substructure velocity distributions are close (< 
300 km s _1 ), then the DS Test cannot always identify real 
substructure. For groups with fewer than ~ 20 members, 
the peaks must be at least 900 km s _1 apart in order for the 
false negative rates to fall below 5 per cent. 

We also find that not only is the DS Test sensitive to 
the e re dahift parameter, but that the level of sensitivity is 
dependent on the number of members in the host group. 
This is best seen by looking at the false negative rates listed 
in the second set of values listed in Table IA4I where we 
use a constant value of Chast = 500 km s _1 for all values of 
"-members- From this section of the table we see that when the 
velocity distribution of the substructure is located at lchost > 
the false negative rates are 43, 49, 34 and 4 per cent for 
mock groups with 10, 15, 20 and 50 members. At 2crhost the 
rates are 10 per cent for groups with 10 and 15 members and 
per cent for groups with 20 and 50 members. This results 
indicates that for groups with fewer member galaxies, the 
velocity distributions of the host and substructure must be 
very distinct in order for the DS Test to detect substructure. 
Alternatively, groups with more than 20 members can have 
substructure galaxies with a velocity distribution embedded 
within the host distribution and still be identified by the 
test. 

A5 Number of Members in the Input 
Substructure (n su b) 

In Table lA"3l we list the false negative rates for mock groups 
with number values of 10, 15, 20 and 50 as a function of the 
number of members in the input substructure (n su b) (see 
Section 2.2.3 for discussion). 

A6 Velocity Dispersion of the Host Group (ohost) 

In addition to the above tests, we also check to see if chang- 
ing the velocity dispersion, and therefore mass, of the host 
group (<Thost) affects the observed false negative rates. We 
present the results in Table lA"6l where we see that for groups 
with more that 20 members the dispersion of the host group 
does not significantly increase the rate of false negatives. 
However, for systems with fewer members and larger values 
of erhost (> 700 km s _1 ), the DS Test is more likely to miss 
true substructure. Fortunately, observed groups with 10 or 
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Table Al. False Negative Rates: Dependency on the Angular Size of the Input Substructure (c p0 sition) 



^position 
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"■members 20 
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48 


26 


20 
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0.5 


83 


74 


73 


27 



"These mock groups have the following input parameters; e pos ition = 0.5 Mpc, o- rcdsll ;f t = 100 km s _1 , e r edshift = 1300 km s _1 and 

°host values listed in Table fT] Only the o- pos j t ion parameter is varied for these trials. 
''Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 mc 

shuffles. 

Table A2. False Negative Rates: Dependency on the Location of the Input Substructure in Position Space (e pos ition) 
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a These mock groups have the following input parameters; <Tp OS ition — 0.01 Mpc, c rc dshift — 100 km s , e rc dshift — 

1300 km s _1 and 

""host values listed in Table [T] Only the e pos ition parameter is varied for these trials. 
^Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 mc 

shuffles. 

so members do not generally have such high dispersion val- 
ues. 
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Table A3. False Negative Rates: Dependency on the Velocity Dispersion of the Input Substructure (o- re[ jshift) 
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"These mock groups have the following input parameters; o" pos ition = 0-01 Mpc, e pos ition = 0.5 Mpc, e r cdshift = 1300 km s — 1 and Oh os t 
values listed in Table [T] Only the oVedsMft parameter is varied for these trials. 
^Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 me 

shuffles. 



Table A4. False Negative Rates: Dependency on the Location of the Input Substructure Along the Line-of-Sight (i.e. redshift space) 

( e redshift) ■ 



e rcdshift ^members 10 ^members 15 ^members 20 ^members 50 
km S _1 
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"These mock groups have the following input parameters; <x pos ; t i on = 0.01 Mpc, e pos ibion = 0.5 Mpc, eredsMft = 1300 km s — 1 and Ohosb 
values listed in Table [T] Only the CT re< j s hift parameter is varied for these trials. 
^Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 me 

shuffles. 

c These mock groups have the following input parameters; <7 pos ; tion = 0.01 Mpc, e pos ition = 0.5 Mpc, e rc dshift = 1300 km s -1 and a 
constant c^ ost value of 500 km s — 1 for all values of n mcm t,ors- Again, only the o- re( jshift parameter is varied for these trials. 
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Table A5. False Negative Rates: Dependency on the Number of Members in the Input Substructure (n su b) 
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"These mock groups have the following input parameters; imposition = 0.01 Mpc, e pos ition = 0.5 Mpc, CT re< j s iuft = 100 km s , 
e redshift = 1300 km s _1 and Oh ost values given in Table [T] For these simulations we only change the number of members in the 

substructure (n su t,). 

^Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 mc 
shuffles. A null entry indicates that no trials were run with the associated n su t, value. 

Table A6. False Negative Rates: Dependency on the Velocity Dispersion of the Host Group (<Th os t) ■ 
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"These mock groups have the following input parameters; CT pos ition = 0.01 Mpc, e pos ition = 0.5 Mpc, o" ro( jshift = 100 km s _1 and 
e redshift = 1300 km s . For these simulations we keep the values of the input substructure constant, but change the velocity dispersion 

of the host group (<T host ). 

6 Values quoted are the false negative rates, given as a percentage, obtained for 100 trials with each set of inputs, using 100 000 mc 

shuffles. 



